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We describe a class of unified theories of electromagnetism and gravity. The 
Lagrangian is of the BF type, with a potential for the B-field, the gauge group 
is U(2) (complexified). Given a choice of the potential function the theory is a 
deformation of (complex) general relativity and electromagnetism, and describes just 
two propagating polarisations of the graviton and two of the photon. When gravity is 
switched off the theory becomes the usual non-linear electrodynamics with a general 
structure function. The Einstein-Maxwell theory can be recovered by sending some 
of the parameters of the defining potential to zero, but for any generic choice of the 
potential the theory is indistinguishable from Einstein-Maxwell at low energies. A 
real theory is obtained by imposing suitable reality conditions. We also study the 
spherically-symmetric solution and show how the usual Reissner-Nordstrom solution 
is recovered. 
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I. INTRODUCTION 



In the seventies Plebanski proposed 



i] 



a reformulation of general relativity where the 
basic dynamical object was taken to be a collection of two-forms satisfying certain algebraic 
constraints. Essentially the same formulation was rediscovered a decade later via a com- 
pletely different path (canonical transformation on the phase space of GR) by Ashtekar . 
A relation between the new Hamiltonian formulation |2[ and Plebanski's theory was eluci- 



dated in 



3j and works that followed, in particular [4]. These developments have made it 
clear that general relativity can be formulated as a theory whose phase space is that of pairs 
(connection, canonically conjugate electric field), i.e. as a gauge theory. This immediately 
suggested that it might be possible to unify gravity with other interactions by considering 
a version of Plebanski- Ashtekar theory written for a gauge group larger than (complexified) 
SU(2) that gives GR. This unification program has been pursued in jf|, |z] at the level 
of the Hamiltonian formulation and in 8] at the covariant level. The almost unknown work 
js] is particularly impressive for it shows how a Plebanski type Lagrangian for the group 
Uc(2) = GL(2, C) describes (complexified) unified Einstein-Maxwell theory. The real theory 
can then be extracted by imposing suitable reality conditions. 

A few years ago one of us has proposed [9[ a class of gravity theories that can be thought 
of as deformations of (complexified) GR in that the new theories continue to describe, as 
general relativity, just two propagating DOF. These gravity theories are most naturally 
formulated using a Plebanski-type Lagrangian where the constraint term for the two-form 
field is replaced by a certain potential term, see below. The idea of these deformations 
is one familiar from e.g. the relation between the scalar field (ft 1 theory with a Mexican 
hat potential and the non-linear sigma model. Here, instead of constraining 7 to lie on 
the sphere ((p 1 ) 2 = const, which gives the non-linear sigma-model, one can allow the ((ft 1 ) 2 
component of the field to fluctuate but give it a large mass. The original theory with the 
constraint ((p 1 ) 2 = const is then recovered in the limit when this mass is sent to infinity. 
Similarly, in the context of theories [9] the fields that were set to zero by the constraints of 
Plebanski theory are allowed to fluctuate but these fluctuations are weighed by a potential. 
The original theory corresponding to GR is recovered in the limit of an infinitely steep 
potential (or at low energies). However, in the gravitational case there is one fundamental 
difference from the scalar field/non-linear sigma-model relationship: replacing the Plebanski 
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constraint term by a potential does not introduce any new propagating degrees of freedom. 
This occurs because the fields that are allowed to fluctuate do not have kinetic terms and 
are thus non-dynamical, see e.g. 10[ for a further discussion of this point. 

As it often happens, the class of theories (9] has been arrived at without knowing that it 
has been considered in the literature previously. Thus, the same class of theories (but in a 
rer different "pure connection" formulation) was described and studied a decade earlier 



rat 
in 



ll| and following works by Bengtsson and Peldan. 

The unification proposal in the context of Plebanski-type theories {9] has been studied 



in 12[ and more recently in |13f j . These references studied the unification that enlarges the 
gauge group of gravity formulated as a real group SO(l, 3) Plebanski-type theory J4]. How- 
ever, in this case, once the Plebanski constraint term is replaced by a potential, additional 
propagating modes get introduced, see [yj and a more recent explicit description in 15|. 
Thus, if one is to avoid extra propagating modes (with some of them with wrong sign kinetic 
terms), there is no choice but to continue working in the chiral (and complex) context of 
the original Plebanski theory fl]. The unification proposal has been studied in this setting 

16). 



in a recent paper by two of us 



The aim of the present paper is to continue the study of unified Gravity- Yang-Mills 
theories of BF type. Here we consider a particularly simple case when the underlying gauge 
group is Uc(2) = GL(2, C). We obtain a simple unified gravi-electro-magnetic theory, and 
our aim is to shed some light on its properties. To this end we first look at the pure 
electromagnetic sector of the model. This is obtained when the gravitational interactions 
are switched off by setting the gravitational fields to their Minkowski spacetime values. The 
resulting theory turns out to be just the most general non-linear electrodynamics with the 
Lagrangian being an arbitrary function of two invariants E 2 — B 2 ,EB. Such models have 
been studied in the literature in the past, see in particular 19) and works by Plebanski and 
co-authors, including [20I ] . The usual Maxwell Lagrangian can be obtained in a limit when 
some parameters of the defining function are sent to zero. 

This behaviour of the theory in the purely electromagnetic setting suggests a useful 
analogy for what happens in the case of deformations of pure GR [9J . Thus, the deformation 
of Einstein's GR to a more general class of theories [9] can be seen to be analogous to the 
deformation of Maxwell's theory to a general class of non-linear electrodynamics theories. 
Indeed, in the Maxwell case the deformation to a non-linear electrodynamics theory consists 
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in adding to the Lagrangian all possible invariants constructed from the field strength. 
Similarly, the gravitational Lagrangian that corresponds to deformations of GR {9] is given by 
a rather general expansion in invariants constructed from the curvature, see 2l|. The arising 
gravitational Lagrangians are not completely general though, as the number of propagating 
DOF described by them is unchanged as compared to GR, a non-trivial property that would 
not hold for a random Lagrangian constructed from curvature invariants. 

Returning to our unified gravi-electro-magnetic theory, the general count of the degrees 
of freedom given in [16J establishes that it is a deformation of both Einstein and Maxwell 
theory with the key property that the number of propagating DOF described by this model 
is unchanged as compared to Einstein-Maxwell. The deformation is controlled by a certain 
potential function, see below, and if one so wishes can be switched off in a continuous fashion. 
Moreover, as we shall explain below, the deformation is only of significance at Planckian 
energies, while for low energies the theory with any generic choice of the defining potential 
is indistinguishable from Einstein-Maxwell. 

We will pay particular attention to the issue of how the real physical theory is extracted 
from an originally complex formulation, the issue that was only touched upon in [if]]. Our 
present analysis of the GL(2,C) model suggests that the prescription for dealing with re- 
ality conditions that was advocated in 16] must be changed. In particular, the somewhat 
unnatural step of taking the real part of the action that was used in [16| in an essential way 
is now eliminated altogether. This implies that the analysis of the Yang-Mills sector given 
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161 ] needs to be changed. This will be presented elsewhere. 



Another possible way to think about the theories that we consider in this paper is that 
they arise by replacing the constraint term of the theory studied by Robinson js| by a 
potential term. Our goal is to better understand the physics of this type of unified theories 
by considering the simplest setting with a clear physical interpretation. 

Our final remarks are as follows. The natural question to confront this work with is: 
Why should one be interested in any deformations of Einstein-Maxwell theory? Indeed, 
this is a well-tested theory correctly describing our world in a vast range of scales. The 
first motivation is aesthetic: The two theories are quite different — one is about spacetime 
metric, the other is about a U(l) connection — and if there exists a reformulation that puts 
them within the same framework (at the price of allowing both theories to get deformed at 
very high energies) it is certainly worth studying. 
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The second motivation has to do with non-renormalizability of GR. The latter implies 
that GR cannot be the gravity theory relevant at Planck energies. The hope is then 9j that 
by sufficiently enlarging the class of gravitational theories one can obtain a class containing 
the sought ultra-violet completion of Einstein's theory. We emphasise that this differs from 
the usual expectation existing in the particle physics community, which is that some new 
DOF (possibly an infinite number of them) become relevant at Planck energies, and resolve 
the high energy problems of Einstein's theory. In contrast, the scenario envisaged by our 
approach is that it is simply the dynamics that changes as one goes to the UV, while the 
dynamical content of the theory (i.e. the number of propagating modes) remains unchanged. 
This can be seen to be a variant of the asymptotic safety scenario 13] for quantum gravity. 

In the case of a unified Einstein-Maxwell theory the motivation is the same. Thus, similar 
to the gravitational sector scenario, the idea is that it is only the dynamics of the Maxwell 
part of the theory that is to change in the UV, while the dynamical content remains intact. 
This is exactly what is achieved by our unified model — our deformations of the Einstein- 
Maxwell theory leave the propagating mode content of the theory intact, only the dynamics is 
changed (in the UV region). The hope is then that this change in the dynamics may capture 
the non-trivial new physics occurring at high energies. We note that this motivation is 
similar to the original motivation behind the Born-Infeld non-linear electrodynamics 18]. 

The organization of the paper is as follows. The next section describes the class of theories 
that we are going to study in this paper. In Section [TTT1 we switch off the gravitational sector 
and study the resulting non-linear electrodynamics theory. In Section |IV] we switch the 
gravitational force back on and study the spherically-symmetric solution of the theory. We 
conclude with a discussion. 



II. THE CLASS OF THEORIES 

In this section we briefly desribe the class of theories that we are going to study for 
an arbitrary complex Lie group G, and then specialize to the case G = GL(2, C) at hand. 
A more detailed discussion of the general G case can be found in jig . One first studies 
a complex theory on a 4-dimensional complexified manifold, with a Lagrangian depending 
holomorphically on all the fields. At a later stage reality conditions are imposed to extract 
a physical real section. 
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A. The class of theories for a general gauge group 

The main dynamical field of the theory is a Lie- algebra valued two-form B 1 , where the 
index J is a Lie algebra one. The other field is a connection one-form A 1 . To write down 
the action, one needs to choose an invariant bilinear form on the Lie algebra of G. When G 
is not simple, this is not unique, and we assume that a choice has been made. Let us denote 
this form by gu. The action of the theory is then of the so-called BF type, with an extra 
potential term for the 5-field: 

S[B, A] = iJ guB 1 AF J - \v{B l A B J ). (1) 

Here i = y— 1 is a factor introduced for future convenience and F 1 is the curvature two- 
form of the connection A 1 . The potential term is non-standard and needs explanation. The 
potential V is a G-invariant, holomorphic, and homogeneous function of order one in its 
arguments. Thus, since its argument is a matrix-valued 4-form X IJ = B 1 A B J , we have the 
following properties of V : V(gXg T ) = V(X),Wg G G; V(aX) = aV(X). The homogeneity 
property is important, for it makes the potential term a scalar- valued 4-form on the manifold, 
which can be integrated to obtain the action. 

A somewhat more practical, but less compact form of writing the action is to introduce 
the density weight one matrix: 

h" := \e^Bl v B J pc . (2) 

Then 

ViB 1 A B J ) = —V(h IJ ) d A x, (3) 

where the homogeneity property makes the right-hand-side to be of the right density weight 
to be integrated over the manifold. We use conventions: dx^ 1 A dx v A dx p A dx a = —W" vpa d^x, 
which explains the minus sign in ([3]). The potential function is then that of ratios of appro- 
priate powers of all the invariants that one can construct from h IJ , see below. 

The field equations that follow from ([1]) are as follows. Varying the action with respect 
to the connection one gets: 



DaB 1 = 0, 



(4) 
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where Da is the covariant derivative with respect to A. This equation can often be in- 
terpreted as an equation allowing one to find the connection in terms of derivatives of the 
two-form field. Varying the action with respect to B 1 one gets: 

1 C^V^ 

9ijFJ = 2d&> (5) 
which, once A 1 is expressed in terms of derivatives of B 1 , can be interpreted as a second- 
order differential equation for the two-form field. To give meaning to the right-hand-side of 
this equation we use the parametrisation of the potential d3J). The field equation becomes: 

9ljF t» = ^n B ^- ^ 

The derivative on the right-hand-side is now the usual derivative of a function of a matrix 
and can be computed without any difficulties. 

B. The case of G = GL(2, C) 

Let us now specialize to the case of interest G = GL(2,C). The Lie algebra in this case 
is 4 (complex) dimensional and splits $j = sc»c(3) ©Uc(l). Up to rescalings, there is a unique 
invariant bilinear form in each factor. Thus, if we split / = (i, 4), where % = 1,2,3 is a 
soc(3) = su<c(2) Lie algebra index, then the most general bilinear form is: 

(X, Y) = k^X'Y 1 + k 2 X a Y 4 , (7) 

where X 1 , Y l , X A , Y A are components of X 1 , Y 1 and 5^ is the usual invariant form on soq(3). 
The curvature components are: 

F l = dA i + ^e) k A j A A k , F 4 = dA\ (8) 

where e*- fc are the soc(3) structure constants. The first BF term of the action then takes the 
following form: 

ik x J 5 ij B i A F j + \k 2 J B 4 A dA 4 . (9) 

Since the normalizations of the two-form fields are not yet fixed we can freely absorb the 
constants Kip into the fields, and we shall do so. 

Let us now discuss the potential term. Let us introduce the following quantities: 

h ij ■= l^B^B^, ft := \~e^B; u B 4 pa , $ := \~e^B 4 u B 4 a . (10) 
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The matrix h IJ is then: 

(< u *|. (id 

The G-invariants of this matrix are: 

Tr(h), Tr(h 2 ), Tr(h 3 ), (0) 2 , & (12) 

where the traces of powers of the matrix h l i are computed using the invariant metric dij, and 
(4>) 2 = Sij(j) l (j)^ . We can take any of these quantities as the basic one, and construct ratios of 
the other quantities and powers of the basic one to form quantities invariant under rescalings 
of B 1 . It is convenient to choose as the basic quantity Tr(/i). The potential function can 
then be written as: 

V(h IJ ) - ^ f ( Tl{h2) Tl{ ~ h3) ^ f t ) (13) 
3 ^(Tr^)) 2 ' (Tr(/i)) 3 ' (Tr(/i)) 2 ' Tr(h) J ' 

where / is now an arbitrary function of its 4 arguments. 

The full action, written in terms of components of forms, is: 

- i J d'x Qe^(%4^ CT + B^F%) - \v{h")j . (14) 

Varying this with respect to the two-form field components one can easily obtain the field 
equations. It is most compact to write them using the form notations: 

F i = ™-B> + ^B\ (15) 
dW <9(0) 2 

dV ~. ■ dV A 
dA = -^ftB* + -LB\ (16) 
<9(0) 2 

where all partial derivatives of the potential can be obtained in an elementary way from 
( fl~3|) . We note that it might appear that a factor of two is missing from the last term on 
the right- hand- side of the first equation, and the first term of the second. However, let us 
carefully compute the variation. We have, dropping unessential constant factors and the 
integral sign: 

^■""(SB^ + SB^) = J|-a« + + gtf . (17) 
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Now, computing the variations on the right-hans-side from the definitions ({TO]) we have: 

f h ij = \i^ 5B %B% 
W = -/ vpa (SB^B^ + B^5B 4 a ) , 
# := \e^°8B%B 4 p(T . 

We now substitute these into ( JT71) and equate to zero the coefficients in front of independent 
variations 5B l pv1 5B pu . We get precisely (TO)]) . ( TTB|) . 

The equations obtained by varying the action with respect to the connection components 
are: 

dB i + e) k A j A B k = 0, dB 4 = 0. (18) 

The first equation here can be solved for the components of A 1 in terms of the derivatives 
of B % . One then substitutes the solution into (1151) and obtains a second-order differential 
equation for B % involving also B . The latter is found by integrating dB 4 = 0, and then the 
connection A 4 is found from (fT6l) . Below we shall see how this procedure works explicitly 
by working out the spherically-symmetric solution of our theory. 

We also note that the equations of our theory are very similar to those of the unified 
theory , with the main difference being that the constraints B l A B^ ~ 8^ and B l A B 4 = 
of are absent in our case. Related to this is the absence on the right-hand-side of 
the Lagrange multipliers that imposed those constraints. Their role is now played by the 
derivatives of the potential function. This is precisely analogous to what happens in the 
case of deformations of pure gravity, where the constraint term in the action is replaced by 
a potential term, and the Lagrange multipliers on the right-hand-side of field equations for 
B l get replaced by dV/dh^ . Thus, the theory that we are considering is a deformation of 
the Einstein-Maxwell theory of precisely the same type as the SL(2, C)-based theory with a 
potential is a deformation of Einstein's GR. Similarly to the case of pure gravity, we shall see 
that it is possible to send some of the parameters of the potential to infinity to recover the 
usual Einstein-Maxwell theory. To understand how this happens, it is useful to first switch 
off the gravitational force, and consider what the theory under consideration becomes as a 
purely electromagnetic theory. 
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III. NON-LINEAR ELECTRODYNAMICS 

A. A version of non-linear electrodynamics 

In this section we switch off the gravitational part of the theory by fixing the soc(3) = 
5u<c(2) part of the 2-form field to be given by: 

B i = T} = idt A dx l - ^e l jk dx j A dx k , (19) 

which corresponds to the Minkowski spacetime background. We further expand the B 4 field 
into the basis of self- and anti-self-dual two-forms: 

B 4 = + = (¥ + ¥)idt A dx l + (¥ - &)-e iik dx j A dx k , (20) 

where and are complex functions, and S l are anti-self-dual two-forms 

t? = idt A dx l + ^e l jk dx j A dx k . (21) 

We now compute the action (fH|) on this field configuration. Using X* A £ J = —1\b % id 4 x 
we get: 

S[$, A 4 } = J d 4 x - WE'^dfrA^ + /($ 2 , $ 2 - ^ 2 )) , (22) 

where / is an arbitrary function of its two arguments. The action depends on fields $, A 4 
that are at this stage all complex. In anticipation of the reality conditions to be imposed on 
the connection A 4 , let us rewrite the Lagrangian in terms of a new connection A: 

A 4 = iA . (23) 

We will later require this connection to be real, with the original A 4 thus being an u(l) 
connection. Using the explicit form of ffT9]) . ( 12 ip we have: 

S[$, tt, A] = J d 4 x ({¥ - ¥){d Ai - d t A ) - i{¥ + &)ej h d s A k + /($ 2 , $ 2 - ^ 2 ))(24) 

It is now clear that the combination ^ — $ l plays the role of the momentum conjugate to 
the spatial projection of the connection A: 

E* := (25) 

and the combination 

Q* : = i(^ + $ l ) (26) 
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is non-dynamical, to be eliminated via its field equation. The action in the Hamiltonian 
form thus becomes: 

S[E, Q, A] = J d 4 x (E%Ai + A d t E* - Q'B, + /((E 2 - Q 2 )/4 + (i/2)EQ, iEQ)) ,(27) 

where we have introduced the magnetic field: 

Bi := eJ h djA k . (28) 

Once the field Q l is eliminated by solving its field equation, we get the non-linear electro- 
dynamics action in the Hamiltonian form: 

S[E, A] = J d 4 x (E*(9 Aj + Aq^E* - H(E, iB)) , (29) 

where H is the Legendre transform of the original potential function / with respect to 
the Q variable. Below we will see how this procedure works explicitly by working out the 
Lagrangian for the function / expanded in powers of its arguments. With the Hamiltonian 
H being a Legendre transform of an arbitrary Lorentz-invariant function, this is the most 



general non-linear electrodynamics Lagrangian, see e.g. [191 ]. [20]. The only difference with 
the Lagrangians typically considered in the literature is that in our case the dependence on 
the invariant EB is with a factor of i = in front, and so it is in general complex even 
after the reality condition A, E e R is imposed. The presence of this extra imaginary unit 
in the action makes the action invariant under a simultaneous operation of parity inversion 
and complex conjugation, similar to what happens in the case of the pure gravitational 
modified theory, see [2l|. This is a very interesting feature of the class of theories considered, 
whose interpretation is still to be understood. In contrast, the non-linear electrodynamics 
real Hamiltonians containing odd powers of EB are, in general, parity violating (if the 
coefficients in front of these terms are taken to be usual scalars). It is however clear that the 
same constraint that is imposed on the Hamiltonian of the usual non-linear electrodynamics 
to have a parity-even theory in our case will produce a real Lagrangian. This will be our 
strategy for dealing with reality conditions below. We leave the more interesting case of 
non-Hermitian Hamiltonians containing odd powers of EB (and its physical interpretation) 
to further research. 

Now, to get a better insight into this theory let us consider its linearisation, in which 
only terms quadratic in the fields are kept. 
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B. Linearised theory 

Unlike considerations of the previous subsection where we have derived the action in the 
Hamiltonian form and kept E* as an independent field, we will now integrate out all fields 
apart from the connection and produce a more familiar Lagrangian that depends only on 
the field strength. At the linearized level we should only keep the terms: 

/(2)($2 ; $2 _ ^2) = + 2 ($ 2 _ ^ (30) 

in the expansion of the function / in Taylor series, where a and 7 are constant parameters. 
Once this is done, we can integrate out the fields $ l , from the action. The solutions for 
$*, are given by: 

& = —E^d^At, ¥ = -Iz^d^Ai , (31) 

a + 7 J 7 J 

and the resulting action is: 

S[A A ] = -\J d'x ^-L-(p^d^At) 3 - ^^At) 2 ^ . (32) 
Using the identities 

jyvY*p* = 2rf^rf^ v - ie^** 7 , f} pv Y} p(y = 2rf- p T]^ v + ie pupa (33) 
that can be checked by an elementary computation we get: 

S M = \ g " ^) / "** + 1 g + ^) / • (34) 

where = c^A^ — d u A A . Thus, modulo the (purely imaginary) second term that is a total 
derivative, we get the following action: 

^i^W^ 1 ^- < 35 > 

Let us note that very little in the above analysis depends on the fact that the gravita- 
tional part of the two- form field was chosen to be (fT9l) . One can see that the procedure 
of integrating out the B A two-form field can be carried out in the same way whenever 
B % A Bi ~ 5^. Thus, whenever the gravitational background is chosen to be "metric", in 
the sense that the Plebanski constraint B % A B^ ~ 5^ is satisfied, it can be seen that the 
linearised electromagnetic Lagrangian is just the Maxwell one, with the metric being the one 
defined by declaring the two-forms B % to span the space of self-dual two forms. This means 
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that the linearised electromagnetic theory is the usual Maxwell electrodynamics not only 
when considered around the Minkowski spacetime, but for any fixed metric background. On 
the other hand, when the condition B % A ~ 8^ is not satisfied (non-metric case using 
the terminology of [9|]), the linearised electromagnetic Lagrangian is different from that of 
Maxwell theory. This means that on a non-metric background light no longer has to follow 
geodesies of the metric defined by B l . Of course, such non-metric backgrounds are only 
of significance in the high-energy regime (small distances). So, we can safely ignore them 
for low energies. Still, it would be interesting to study the effects of non-metricity on light 
propagation; we leave this to further research. 



C. Linearised reality conditions 

Assuming (for simplicity) that both a, 7 are real and positive, we easily deduce the 
linearized level reality conditions that must be imposed on our fields. Thus, the condition 
that A 4 is purely imaginary, which is appropriate if we want to think of A 4 as the u(l) 
component of a connection field, gives the correct Lorentzian signature action. Thus, for 

A 4 = iA, AgR, (36) 

we get: 



5[A] = -_- / d 4 xF^F^, (37) 

4 #u(l) J 

where = d^A^ — <9„A M is the fields strength and the coupling constant is: 

9l(i) = T(« + 7)/«- ( 38 ) 

We now note that, if desired, we can obtain the [8| version of the electrodynamics in 
which the field B 4 is purely anti-self-dual by sending a — » 00. Indeed, as clear from (I3T]) . in 
this limit that describes the self-dual part of B 4 goes to zero. In this limit the coupling 



constant of our Maxwell theory becomes g^,-^ = 7. Thus, the theory considered in is 
easily recovered. 

When a —7- 00 the two-form field B 4 becomes (proportional to) the anti-self-dual part of 

the real field strength F^ v . In the case of finite a the reality conditions that B 4 satisfies are 

much more involved. We get: 

4 _ a + 27 F a i P a F m 

^ 7 ( a + 7 )^ + 7 ( a + 7 ) 2 e ^ tp °- [6J) 
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The structure arising is typical for the theories under consideration in that the part of the 
expression that carries the e^ vpa tensor contains an additional factor of i as compared to the 
part that does not contain e pvpa . 



D. Non-linear electrodynamics 

Above we have analyzed the theory with the potential function / truncated to its 
quadratic terms in the quantities $ 2 , \l/ 2 . To understand the structure of the full non-linear 
theory we expand the potential and keep higher powers of $ 2 ,\l/ 2 . Thus, let us see what 
happens at the next order, which is quartic (Lorentz invariance prevents us from having any 
cubic terms). The quartic order part of the potential can be parametrized as: 

/ (4)($2 ) $ 2 _ = X ($2) 2 + ^2^,2 + !^ 2)2; (4Q) 

where Xi & an d £ are constant parameters. 

One can now vary the action fl22|) with respect to $ l , ty l and solve for these fields perturba- 
tively in powers of A 4 . We get for the cubic order terms: 

=j^^At (j^y^dAr + ^dAr) , (41) 

V {3)1 = - ^p^d^Al ( 6 . 2 (£M 4 ) 2 + ^(ZdA 4 ) 2 ) , 
T V(« + 7) 7 / 

where we have introduced a compact notation 

(EcM 4 ) 2 := ^dX^^K = \^ V K ~ -/ VP °FU F U ( 42 ) 

and similarly for (£cL4 4 ) 2 . Now we can compute the quartic order Lagrangian, with the 
result being: 

£<4> = ^if^ff + 2l Ha + 1 / WAimdAi)2 + 4^ ((SM4)2)2 - (43) 
This can be expanded in terms of the usual field strength invariants. Thus, using 

(e^F* u Fp 2 = -8(F^Fp 2 + lQF^F^F^F^ , (44) 
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we get the following Lagrangian: 



16 1 ^ V(« + 7) 4 7 2 (« + 7) 2 7V 1 J 



.^p^F^F^F^ I x 2S < * 



4 " " p CT V(« + 7) 4 7 2 (« + 7) 2 1 A 



We can now substitute here the linearized reality conditions (|36|) and obtain the Lagrangian 
for the real-valued connection. However, we note that now, unlike what happened in the 
quadratic order of the theory, the imaginary term in the Lagrangian is no longer a total 
derivative. Thus, as we have already discussed above, the non-linear action for the real 
connection (J36J) is, in general, complex. This is precisely similar to what happens in the case 



of the effective gravitational Lagrangian, see 



21] . There the metric Lagrangian that one gets 



from a similar BF-type theory with a potential (but in the case of G — SL(2,C)) at cubic 
order in the curvature in general contains an imaginary term that is not a total derivative. 
Similar to what we are seeing here, in the purely gravitational case it is also the higher- 
order interaction term that is in general complex, while the theory linearized around the 
Minkowski background does not exhibit any complexity issues. We also note that, similar to 



what happens in the case 



2l| of pure gravity, the imaginary term is odd under parity. Thus, 



the full Lagrangian is invariant under the operation of complex conjugation accompanied by 
parity inversion. 



E. Reality conditions 

There are several strategies that one could follow when facing such a non-Hermitian 



Lagrangian. One, advocated in jl6j, [2l| is to impose the linearized reality conditions and 
take the real part of the full non-linear action. As was, however, realized more recently in 
the context of work [22j on the purely gravitational theory linearized around the expanding 
FRW background of relevance for cosmology, this real part of the action prescription does 
not in general produce a consistent theory. In the case of the FRW background the problem 
arises when one considers the gravitational waves (tensor perturbations). 

Let us describe what the problem is in some more details. As in our electromagnetic 
considerations above, in the purely gravitational case one starts from the complex action 
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and "integrates out" the non-dynamical fields to obtain an action that depends only on 
the physical metric. The action one gets is a functional depending holomorphically on 
the complex "physical" fields. One has to impose some reality conditions to extract the 
real action. Since the action depends on the complex fields holomorphically, one can instead 
consider the, say, real part of the action, and vary it with respect to real and imaginary parts 
of the fields. The arising Euler-Lagrange equations are the same as the real and imaginary 
parts of the complex Euler-Lagrange equations one obtains from the holomorphic action (this 
follows from Cauchy-Riemann equations). Thus, the real part of the holomorphic action 
considered as a functional of real and imaginary parts of all the fields carries exactly the 
same information as the original holomorphic action and can be taken as the action for the 
theory. However, this action depends on twice the number of physical fields, and, moreover, 
kinetic terms for the "imaginary" parts of the fields are typically negative-definite. Thus, this 
action describes twice the number of propagating modes of the physical theory, and is badly 
unstable. Reality conditions are needed to select a good physical sector of the theory, which 
describes half the modes of the complex sector and is void of any instability problems. It is 
natural to require that the reality conditions one imposes are some second-class constraints 
that cut the dimension of the phase space by half. However, as a consideration of simple 
examples shows, for an action that is obtained as the real part of a holomorphic action, it is 
in general not consistent to impose the constraint that the field is real. The reason for this is 
that the condition that this constraint is preserved in time generates a secondary constraint, 
and the condition that the secondary constraint is preserved in time in general produces a 
new constraint that is not equivalent to the original constraint of the reality of the field. 
Thus, in general, requiring the field to be real imposes more constraints than one would 
want. So, in general the dynamics of the Lagrangian such as f )45|) (or the real part of this 
Lagrangian with all fields complexified) is not consistent with the reality condition that the 
physical field is real. If one imposes this condition at some instant of time (and arranges the 
momentum to be real as well), the dynamics will in general generate an imaginary part of the 
field. So, the strategy of dealing with the problem of reality conditions for non-Hermitian 
Lagrangians should be more sophisticated. We will leave any attempt at such to further 
research. 

For the purposes of this paper we note that, at least in our case of non-linear electro- 
dynamics, we can restrict the potential function defining the theory so that the arising 
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Lagrangian is real (for real connections). This is precisely what is usually done in the con- 
text of non-linear electrodynamics theories studied in the literature, where there is typically 
a restriction on the class of defining functions so that the theory is parity- invariant. Thus, 
in the case of our Lagrangian ( )45j) we can arrange the coefficients in the expansion of the 
function /($ 2 , $ 2 — in such a way that the coefficient in the last term in (H5i) is zero. We 
can arrange things so that no imaginary terms arise in the higher orders of the expansion 
either. This will produce a consistent theory with a real Lagrangian as far as the electro- 
magnetic sector is concerned. Whether a similar prescription is possible in the gravitational 
sector is beyond the scope of this paper. However, in the next section we shall see that at 
least in the spherically-symmetric situation the reality condition that the metric is real is 
completely consistent. 



IV. SPHERICALLY-SYMMETRIC SOLUTION 

In this section we obtain and analyze the spherically-symmetric solution of the gravi- 
electro-magnetic theory described above. As we have already noted above, on non-metric 
backgrounds where B % A ^ 5 l i the coupling of electromagnetism to gravity as prescribed 
by our theory is different that in the Maxwell case. Thus, there are two different ways 
that electromagnetism can be coupled to deformations of GR [9]: one way is to couple 
the electromagnetic potential to the metric defined by declaring the gravitational sector 



two-forms B % to be self-dual. Such a coupling has been studied in [23], where also the 
spherically-symmetric solution was analysed. A different coupling is given by our theory. 
Thus, the spherically-symmetric case field equations that we shall analyse are distinct from 
those in {23 1. We emphasise that this difference is only of relevance for very small scales (or 
Planckian curvatures). For low energies (large distances) the sperically-symmetric solutions 



of both 23] and this work become indistinguishable from Reissner- Nordstrom. 



A. The spherically symmetric ansatz 



We start by making an ansatz for all the fields as dictated by the symmetry. The gravita- 
tional su{2) sector B-fields can be selected as in the purely gravitational case first studied in 
24j . This reference has worked in spinor notations and used a complex null tetrad. However, 
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it is not hard to repeat the analysis for a real tetrad for the usual spherically-symmetric 
metric asatz 

ds 2 = -fdt 2 + g 2 dr 2 + r 2 d6 2 + r 2 sin 2 (#)# 2 , (46) 
where as usual /, g are (real) functions of the radial coordinate r only. The starting point 



of the analysis is to construct the self-dual two-forms for this metric, see e.g. [25] for a 
description of this procedure for the case of Einstein's GR. The modified theory ansatz is 
then obtained by allowing for an extra functions of the radial coordinate multiplying the 
metric B-field ansatz of the GR case. Using the available coordinate freedom one can put 
the B-field in the following convenient form: 

B l = b(ifrdt A d6 - gr sin 9d(f) A dr), 

B 2 = b(ifr sin 9dt Ac/0 - grdr A d9) , (47) 
B 3 = ifgdt A dr - r 2 sin 9d9 A d(p , 

where b is a function of the radial coordinate. When 6 = 1 one gets the usual metric self- 
dual two- form ansatz of relevance for Einstein's GR, see j^ . As we already mentioned in 
the previous section, when the parameters of the potential of the electromagnetic sector are 
chosen so that the purely electromagnetic Lagrangian is real, the metric also turns out to 
be real. Thus, in the spherically-symmetric case one can assume that the metric is real from 
the start. We therefore assume the functions /, g and b, as well as the coordinate functions 
t, r, 9, (ft to be real. 

The ansatz that we make for the B 4 two-form field is a general combination of the 
"electric" and "magnetic" two-forms: 

B 4 = -2cr 2 sin 9d9 A d<f> + 2i fgm dt A dr , (48) 

where c, m are functions of r only, and the numerical constants are introduced for future 
convenience. No reality conditions on c, m are assumed at this stage. 



B. 5-compatible GL(2, (C)-connection 

We now solve 

DaB 1 = dB 1 + fj K A J A B K = 
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for the connection. The gravitational su(2) part of this "compatibility" equation reads: 



D A B i = dB l + e) k A j A B k = 0, 



(49) 



which gives: 



A 1 = - — sin 6d(j), 
bg 

A 2 =^-d9, 



bg 

if 

9 



(50) 



(bfr)' 1 
bfr b 2 r 



dt + cos 9d(fi . 



For the u(l) part of the compatibility equation we have DaB 4 = dB 4 = 0. This implies: 

(cr 2 )' = 0. (51) 

Note that we cannot solve this equation for A 4 , so we will need to find the electromagnetic 
connection from another equation. For now, we make the following spherically-symmetric 
ansatz for it: 



iadt + ip cos 9d<p , 



(52) 



where a is, at this stage, arbitrary functions of r, and the imaginary unit is introduced in 
the expectation that later the reality condition will be imposed requiring the connection to 
be purely imaginary, as appropriate for a u(l) connection. The spherical symmetry requires 
p to be a constant (proportional to the magnetic charge of our system). 
We can now compute the curvature 



F 1 = dA 1 + -f JK A J A A 



K 



(53) 



of the connection that we found (or made an ansatz for) above. We have for the gravitational 
sector 7 = 1,2, 3: 

1 

bg, 



F bg 
/ = <! ^ 



1J 


fi/ 


I (bfr)' 


1 " 


bg 


u 


bfr 


b 2 r 


ij 


'i/ 


'(bfr)' 


1 " 


bg 


.9 


bfr 


b 2 r 



sin 9dr A defy > , 



dtAd9 + bg 



sin 9dt A d(p - bg [ — ) dr A d9 \ , 



(54) 



(bfr)' 
bfr 



1 

b 2 r 



dt A dr 



b 2 g< 



sin #<i# A d<f> : 



and for the electromagnetic field strength 1 = 4: 

F 4 = dA 4 = -m'dt Adr- ipsm(9)d9 A # . 



(55) 
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C. Field Equations 

The remaining field equations to consider are given in ( |T5l) above. Defining the matrix 
h IJ via: 



B 1 A B J = h IJ (-2ifgr 2 sin 6 dt A dr A d6 A d<f>) 



(56) 



we get: 



JJ 



b 2 
6 2 








1 c + m 



c + m 4cm 



(57) 



We can now compute the derivatives of the potential (fTBl) needed in (fT5|) . We get 



<9V 
dp 



/ - 2(26 4 + 1) 



fx 



3(2b 2 + 1) 



, 2fe 6 + l 2(c + m) 2 



(26 2 + l) 3 JJ 3(2& 2 + l) 5 



4cm 



3(26 2 + 1) 



3(26 2 + r 



(2b 2 + iy- 



3(2b 2 + 1) 



dV 

dip 



fi 
3 



(58) 



Here f' n is the derivative of the function / with respect to n-th argument evaluated at 
= diag(6 2 , b 2 , 1), 4> 2 = (c + m) 2 , ip = 4cm. 

It turns out to be very convenient to separate the trace and the tracefree parts in the 
gravitational part, and introduce a separate notation for the electromagnetic part potential 
first derivatives. Thus, let us write the matrix of the first derivatives of the potential as: 



dV 
dh 11 



A-/3 

A-P 

A + 2(3 a 

a p 



(59) 
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where A, /3, p, a are functions of b, c, m given by: 



a / r 4(& 2 -!) 2 
3 Jl 9(26 2 + 1) 2 



^ 2(6 2 - 1)(6 4 - 1) _ 2(c + m) 2 



3(26 2 + l) 3 



s 3(26 2 + 1; 



/I 



4cm 



3(26 2 + 1)' 



9(26 2 + 1] 
c + m 



: 3(26 2 + l) 2 ' 



fi 



3 3(26 2 + l)' 

The field equations (115j) then read, in the gravitational sector: 

A-/3. 



1 

bgr \bg 



b 2 g 2 r 



(bfr)' 
bfr 
1 



1 

b 2 r 



f_ 

fg [g 



(bfr)' 
bfr 



b 2 g 2 
1 

b 2 r 



-- A + 2/3 + 2c(T , 
A + 2/3 + 2mcx . 



The electromagnetic sector equation gives the following two equations: 

a = —fg(a + 2pm), ip = r 2 (a + 2pc). 



(60) 



(61) 
(62) 
(63) 

(64) 



Before we analyze these equations let us describe a convenient change of independent func- 
tions that will eventually allow us to integrate the system. 



D. Legendre Transformation 

As was done in the purely gravitational spherically-symmetric case treated in [24], we 
can think of A as the Legendre transform of the function /. In fact, we have 

/ 

A = 3 ~ X ^ ~ W ~ Zp) ^ 



where 



Thus, we get 



1 — b 2 c + m 4cm 

% = 2—— , y = 2—— , z = — -7; . (66) 

26 2 + l' y 26 2 + l' 26 2 + l v ; 



# A 9A * 9A 

A, :=- = -*, A a := — = -y, A p := - = -z. (67) 

We can use these relations to express the original functions b, c, m appearing in our two-form 
field ansatz in terms of derivatives of the new function A = A(/3, a, p). We get: 

2 2 + Ap c + m = 3A - 2cm = ^ (68) 

2(1 -A,)' C + m 2(1 -A,)' ICm 2(1 -A,)" ^ 
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This gives, for c and m: 

C =-4^ + \/(r3A^ + ^^)J. (69) 

m =-? f A - ,/ A ° | 4A " ) 
i\l-Af, V (1 - A «) 2 3(1 - A »)/ ' 

We have chosen the solution such that m = for A p = 0. Thus, one can now take the 

viewpoint that the theory is parametrized by the function A = A(/3, er, p), and that the 

above relations give us the functions b, c, m once the quantities (3, a, p are solved for. This 

change of viewpoint will allow us to integrate the field equations. 



E. Bianchi identities 



A very powerful method for analyzing the system of equations that we have obtained is by 
rewriting them as differential equations for the functions (3, a, p. These are nothing but the 
Bianchi identities obtained from the equation DF 1 = 0. Alternatively, these equations can 
be obtained directly from the field equations (!6T!) - (!63l) . Thus, differentiating the equation 
fl62|) . and using one of the equations in flBTl) . as well as fl5Tl) we get: 

6/3 



A' + 2/3' + lea' = -- 



(70) 



Another Bianchi identity is obtained by differentiating the second equation in ( IBTI) . We 
have: 



b 2 g 2 r 



{bfr)' 
bfr 



1 

b 2 r 



bg r J b 2 g 2 r 



(bfr)' 
bfr 



1 

b 2 r 



A' - 13'. 



(71) 



We now rewrite the equation fl63|) expanding the terms on the left and dividing the whole 
equation by b 2 r. We have: 



b 2 g 2 r 



(bfr)' 
bfr 



1 

b 2 r 



f 9 



1 



b 2 g 2 r 



(bfr) 
bfr 



1 

b 2 r 



A- (3 + 3(3 + 2ma 
b 2 r 



■ (72) 



Using the second equation in (1611) we express A — (3 in terms of other quantities and then 
take this term to the left-hand-side of the equation. We get: 



1 


I {bfr)' 


1 " 


( l - f ' + s ') 

\b 2 r f gj 


1 


Kb fr)' 


b 2 g 2 r 


bfr 


b 2 r 


b 2 g 2 r 


bfr 



1 

b 2 r 



3(3 + 2ma 



b 2 r 



We now subtract (175)1 from (ITTI) . We get: 



b 2 g 2 r 



(bfr)' 
bfr 



1 

b 2 r 



(b 2 fg)' | b 2 - 1 

b 2 fg b 2 r 



A'-/3' 



3/3 + 2ma 
b 2 r 



(73) 



(74) 
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We should now note that the following equation is true: 

ib 2 fg)' 1 - b 2 



(75) 



b 2 fg b 2 r ' 

Indeed, this is just a rewrite of the first equality in flfJTj) . Therefore, the quantity in the 
second brackets on the left-hand-side of (1711) is zero and we get: 

A' - j3' = W+J™ . (76) 
err 



Equations fl70|) . fl76|) . together with (J5TJ), after the functions b, c, m are expressed in terms 
of (3, a, p via (EH]), become 3 first order differential equations for the 3 unknown functions 
(3, a, p. Once these are found, the electromagnetic connection is found from flMl) . and the 
metric functions g is found from 

' 1 - (A + 2/3 + 2ca)r 2 , (77) 



(bg) 2 

which is easily obtained from (|62|) . The metric function / is then found from (IT5]) . 



F. Consistency 

Yet another Bianchi identity can be obtained from dF A = 0, and is equivalent to the 
statement that the magnetic charge p = const. On the other hand, the second equation in 
( )64|) expresses the magnetic charge in terms of other functions. Differentiating this equation, 
and using the known from (15T!) derivative of c we get: 

a' + 2p'c + 2a /r = 0. (78) 



This equation can be shown to follow from the two Bianchi identities flTOl) . fl76j) and the 
relations (!68|) . Thus, let us show that (17H|) together with (176|) imply (I70I) . We multiply ( I7H|) 
by 2m and, using the last two identities in f )68|) write the result as: 

4ma _ 3p'A p + 3cr'A 



« , , 2a'c. (79) 

r l-Ap 

We now use this, as well as the first identity in f )68|) . to write —2b 2 times f J76|) as: 

2 ± A g+ 3p-A p + 3 < /A, +2 ^ = _6g 
1 - A/3 1 - A/3 r 

The first two terms on the left-hand-side combine to: 

2/^ - 2P'A P + A' - A'Ap 

l-A/3 

Thus, f lHUj) is just flTUj) and the obtained system of equations is consistent 



A' + 2(3'. (81) 
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G. Non-metric gravity 
In the limit A CT = A p = the electromagnetic part of the theory is switched off and 



we recover the spherically symmetric solution 24| of non-metric gravity. The two Bianchi 



identities (170)) . ( 1751) in this case coincide and give the following equation 



(A„ + 2)3' = -^- (82) 



are determined from (177)) , ( I75p . For more 



for 8. After this is solved the metric functions 
details on the pure gravity sector solution see |24j. 



H. Reissner-Nordstrom solution 

Let us now see how the usual Reissner-Nordstrom solution of GR coupled to Maxwell 
can be recovered. First, we should switch off the gravity modifications, which is done by 
putting A^ = which gives b 2 = 1 and the gravitational part of the two-form field becomes 
the usual spherically-symmetric triple of metric two-forms. The simplest way to get the RN 
solution is to set A CT = so that m = —c and the B 4 field ( 14"8)) is anti-self-dual. However, 
let us see the appearance of the charged solution in full generality. We will also allow the 
magnetic charge to be present, to illustrate how the issues of complexity should be dealt 
with. 

First, we need to perform the Legendre transform of the original defining function /. In 
a previous section we have seen that in the absence of gravity modifications, and in the 
case which gives the usual Maxwell theory, this function is given by ( 130)) . Inspection of ( 1571) 
reveals that we have to replace $ 2 — > (c + m) 2 , $ 2 — \l/ 2 — > 4cm. Thus, in the case that gives 
Maxwell theory our defining function is 

f(c,m) = -{c + m) 2 + j{Acm). (83) 

We then easily find a, p from ( 160)) : 

df / \ a(c + m) df 7 

a = WT m J 2{c + m) = ^^^ P = W^) = 2- {) 

The Legendre transform (165)) now gives: 

A = -J?V. (85) 
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Note that this is independent of p, as the original function was linear in p. However, the 
derivative A p cannot be considered to be zero because it must satisfy the last equation in 
( I68p . Thus, in this case the parametrization by A is somewhat degenerate. This can be dealt 
with by declaring the last equation in ( 1681) to be satisfied by definition. This degeneracy is 
removed when one considers more complicated, non-linear dependence on cm. 

We can now proceed to solving the equations. We first find a from the second equation 
in ( jMJ) . Using the value of p given by (|84|) we get: 

(86) 



We now find m from the second equation in fl68|) and get: 

2ip 27 + a 

m = — - - c , (87) 

or a 

where we have used = 0. We now use the first equation in ( 164]) and, in anticipation that 
no modification to the electromagnetic potential will be introduced, put a' = —q/r 2 , where 
q is the usual (real) electric charge. This allows us to express the quantity c in terms of q,p. 
Using fg — 1 (which follows from fl 75 j) ) we get: 

2 7 r 2 ( 7 + a) K J 

Note that this does have the required 1/r 2 dependence on the radial coordinate. The above 
expression for c gives the following expression for a, m 

a= a^ + ip) = (2 7 + a)g - cap 

2r 2 (7 + a) 2 7 r 2 ( 7 + a) v ; 

We note that all quantities a, m, c became complex, so it is by no means obvious that one 
will arrive at a real metric at the end. Note also that, interestingly, the quantities c,m can 
be obtained one from another by exchanging ip <H- q. 

We can now solve for the unknown function (3 using e.g. ( !?0l) . Using A' = — (4/3a)aa' 
and putting all the terms not containing /3 to the right-hand-side we get: 

p, = _ a(q + ip)((2 7 + 3a)q - (4 7 + 3a)ip) 

r 6'jr 5 (a + 7) 2 

The solution with correct behaviour at infinity is: 

a_ r s a(g + ip)((27 + 3a)g- (47 + 3a)ip) 

1 2r 3 6 7 r 4 (a + 7 ) 2 ' 1 ' 
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Note that this quantity is, when p ^ 0, complex even when the reality conditions are 
imposed. 

We can finally find the metric functions from (I77|) . The above analysis does not seem to 
make it plausible that the arising function g can be real. However, once we substitute all 
the quantities we have found above into ( 1771) we obtain: 



Thus, the metric is the usual real Reissner-Nordstrom black hole with electric and magnetic 
charges provided we choose a, 7 so that: 



which is exactly the condition expected from the formula (I38p for the coupling constant. 

To summarize, the analysis of this subsection confirms that there exist a two-parameter 
family of potentials giving rise to unmodified Einstein-Maxwell system. It also illustrates 
how non-trivial can the issue of reality become. Indeed, we have worked with complex 
quantities at intermediate stages of the computation, but at the end all the complexity 
disappeared to give rise to the real metric functions. This could have been expected from 
general considerations, since we have switched off the gravitational and Maxwell sector 
modifications. However, it is reassuring to see this happening explicitly. 

Any departure from the simple choice of the defining potential considered in this subsec- 
tion produces a modified theory, where one can either modify the gravitational sector, or 
electromagnetic, or both. At this stage of the development of the theory it is unclear what 
constitutes a physically interesting potential function. Also, as we have emphasised above, 
our theory only becomes different from Einstein-Maxwell in the region of high curvatures. 
In this regime other effects (quantum mechanics) also become important. Thus, it is not 
at all clear if in this regime any physics can be extracted by looking at the purely classical 
solution. Thus, we refrain from an analysis of the modified black hole solutions at this stage. 
At the same time it is gratifying to know that the theory is simple enough that the problem 
of determining such a solution for a general defining potential reduces to three first order 
ODE's for the functions (3, a, p. 
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-2 




(92) 



a 



(93) 



7(0 + 7) 
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V. DISCUSSION 

In this paper we described a simple class of unified gravi-electro-magnetic theories. The 
Lagrangian of our unified theory is of the BF type, with a potential for the B-field added. 
The gauge group on which the theory is based is (complexified) U(2) = SU(2) x U(l), and 
the SU(2) sector of our model is responsible for the gravitational interactions while the U(l) 
sector describes electromagnetism. 

The Lagrangian of our theory can be argued to be the most general functional of the Lie- 
algebra valued two-form field B and the connection A, subject to the conditions of gauge and 
diffeomorphism invariancdj. As such, it is a very natural starting point for a unified theory 
that describes both gravity and a gauge theory. It is gratifying to learn that the simplest 
possible nontrivial such a theory describes gravity, and the next most simple case describes 
a unified gravi-electro-magnetic theory. Indeed, the simplest possible gauge group that one 
could try to use in the general Lagrangian ([!]) is the Abelian U(l). However, in this case the 
only potential that one could write is B AB, which produces a topological theory. Thus, the 
U(l) gauge group does not lead to a physically interesting theory. The next simplest group 
is SU(2), and for this choice we do get a non-trivial gravitational theory describing as GR 
just two polarizations of the graviton. Increasing the level of complexity one step further 
one takes the gauge group to be U(2) = SU(2) x U(l) studied in this paper. As we have 
seen, this is a unified theory of electromagnetism and gravity. It is quite encouraging that 
the first two simplest non-trivial cases that one encounters when studying diffeomorphism 
invariant gauge theories (OQ) are the two most physically important theories. 

We would like to emphasise that the only assumptions that were used in our choice 
of the Lagrangian are those of diffeomorphism and gauge invariance. It is true that the 
potential function in (pQ) can be arbitrary, so one has not a single theory, but the whole 
class parameterised by the choice of V(-). However, we have also seen that for any choice of 
the potential the theory expanded around the Minkowski spacetime background is for low 
energies indistinguishable from the Einstein-Maxwell theory. Indeed, in the lowest, quadratic 
order in the expansion the gravitational and electromagnetic parts are uncoupled. One can 
then refer to the pure gravity case results in e.g. [l6] to see that the low-energy gravitons 
1 An alternative equivalent viewpoint is obtained by integrating out the two-form field B and considering 
the "pure connection" formulation, in which the theory becomes just the most general diffeomoprhism 
invariant gauge theory. We refrain from considering this point of view in the present paper. 
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are unmodified. On the other hand, the electromagnetic sector is described at low energies 
by the quadratic terms fl30l) which we have seen to lead to the usual Maxwell Lagrangian. 
Thus, for any choice of the potential function defining the theory (apart from that giving 
the topological BF theory) the linearized excitations are the usual gravitons and photons. 
One can also see that the interactions between these are correctly described by noting, as we 
already did above, that on an arbitrary metric background described by two-forms satisfying 
B l A Bi ~ 5^ the linearised electromagnetic sector coupling to gravity is just that of the 
Maxwell theory. Thus, for low electromagnetic field strength and small curvatures in the 
gravitational sector our theory is just the Maxwell theory coupled to gravity. 

Thus, most of the parameters of the potential function V are invisible at low energies, 
and we only have access to the usual Newton's constant and the Maxwell theory coupling 
constant (unmeasurable unless a coupling to matter is introduced). The departure of our 
theory from Einstein-Maxwell would only be visible at high energies (or high field strength). 
Since the only dimensionful parameter in the theory is the Newton's constant, the only 
relevant energy scale is the Planck scale. Thus, the departure of our unified theory from 
Einstein-Maxwell is only going to be visible close to the Planckian domain. For low energies 
that we have access to any choice of V leads to the usual Einstein-Maxwell system. 

We could also have phrased the above discussion in terms of renormalisation group flow 
arguments. Since our theory was seen to describe the degrees of freedom of the Einstein- 
Maxwell theory, and has the same symmetries, we are guaranteed that at low energies 
the theory will be indistinguishable from Einstein-Maxwell, since by general arguments the 
Einstein and Maxwell Lagrangians are the only terms of the effective field theory Lagrangian 
of gravity and electromagnetism that survive at low energies. This is confirmed by our 
analysis. 

The choice of the potential therefore only matters in the high energy regime where the 
theory starts to deviate from its low energy Einstein-Maxwell limit. However, here other 
effects are becoming important as well, of which the most important is of course quantum 
mechanics. Indeed, at high energies we can no longer treat our theory as classical. This in 
particular implies that we are no longer free to specify the potential V, as it will become 
running with the energy scale. So, at high energy where a choice of V would be of importance 
we are no longer free to make this choice, and it is the renormalisation group flow which 
will determine V. This is gratifying, because it means that there is no freedom at all in 
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our choice of the Lagrangian, and it is completely specified by symmetries and the quantum 
behaviour. Unfortunately, we are not able to say anything about the latter at this stage of 
the development of the theory. 

Let us finish with a quick discussion of the open problems of this approach. First, while 
there is some scope for gauge-gravity unification in the context of our diffeomorphism in- 
variant gauge theories, there are other field species in Nature - fermions. These usually 
require a metric (or a tetrad) if they are to couple to gravity, and so it is not at all clear 
that they will be possible to describe in an approach that trades the spacetime metric for a 
collection of two-forms. This is a difficult question that will be dealt with elsewhere. One 
other very important question is that of reality conditions for the whole theory (including 
the gravitational sector). There are some puzzles here associated with the fact that the 
Lagrangians one naturally obtains in our approach are non-Hermitian (this non-Hermiticity 
is only visible in higher-order interaction terms, see e.g. the last term in (l4"5j) . and is of no 
significance at low energies). We will return to this problem elsewhere. For now let us just 
mention the fact that it is now known that the condition of Hermiticity is overly restrictive 
and there are some systems in Nature that are described by non-Hermitian Hamiltonians, 
see 26j for a review. Finally, the most important open problem of the whole approach is 
to study the quantum mechanical behaviour of our theories and show that they continue to 
make sense also as quantum theories. 
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